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Spherically symmetric gravitationally bound, oscillating scalar lumps (boson stars and oscillatons) 
are considered in Einstein's gravity coupled to massive scalar fields in 1 + D dimensional de Sitter- 
type inflationary space-times. We show that due to inflation bosons stars and oscillatons lose mass 
through scalar radiation, but at a rate that is exponentially small when the expansion rate is slow. 
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I. INTRODUCTION 

Current physical field theories explaining the presently 
known fundamental interactions are remarkably success- 
ful to account for phenomena both on particle physics 
(i.e. 10~ 18 m ) and on cosmological (10 26 m) scales, usu- 
ally incorporate various fundamental scalar fields (Higgs 
bosons, axion, supersymmetric partners of fermions, in- 
flaton, waterfall fields etc.). It has been observed since a 
long time ago that self-interacting scalar fields admit at 
least up to D — 4 spatial dimensions spatially localized, 
oscillating lumps, nowadays called oscillonsyj-yj. Oscil- 
lons are reminiscent of spatially localized, time-periodic 
breathers known in one dimensional Sine-Gordon theory. 
There is, however, a fundamental difference, in that os- 
cillons are only approximatively time-periodic since they 
radiate, thereby losing continuously their energy, albeit 
very slowly [||-|g|. Moreover oscillon-type objects exist in 
a broad class of field theories containing in addition to 
scalar, vector and possibly other fields [7, 8]. A promi- 
nent theory admitting classical oscillons is the bosonic 
sector of the Standard Model [9l-fll|. Einstein's grav- 
ity coupled to a free, massive real Klein-Gordon field, 
also possesses oscillon-type solutions, which go under the 
name of oscillatons |12 - tl7J . 

An important problem is to investigate the possible 
pertinence of oscillatons for cosmology where according 
to current theories some scalars fields are necessarily in- 
volved. In recent investigations the evolution of oscillons 
in an inflating background has been investigated |18l - l23j . 
It has been found that an inflating background also in- 
duces a radiative tail in oscillons, distinct from the one 
present in flat space-times, leading to an additional leak- 
age of energy. However this decay rate is also expo- 
nentially suppressed, so inflating oscillons may survive 
for an exponentially long time provided that the horizon 
is far larger than the width of an oscillon. As recently 
shown, due to quantum effects, the decay rate of oscil- 
lons becomes power law-like and in many (but not in all) 
cases dominates over the classical decay rate [24| . There 
are indications that oscillons can be copiously produced 
from thermal initial conditions and a substantial frac- 
tion of the energy of the system will be stored by them 
HM HDI l25[ ■ Therefore it is important to investigate the 
physical role of such long living localized objects. In a 



cosmological setting it is necessary to take into account 
the influence of gravity, i.e. instead of simply considering 
oscillons in a background metric, to investigate solutions 
of the coupled Einstein-scalar system - oscillatons - with 
an expanding, asymptotically homogeneous metric. 

In this paper we undertake a thorough investigation of 
localized, oscillaton-type solutions in Einstein's gravity 
coupled to a massive Klein-Gordon field, when space- 
time undergoes inflation, driven by a negative pressure 
cosmological fluid - i.e. including a cosmological con- 
stant, A > 0. (For recent reviews on inflation theory 
see e.g. [UH.) 

Our results show that once admitting massive scalar 
fields, both in present day cosmology and in the early 
Universe, spatially localized, very long living oscillatons 
exist. The inflaton field generates an energy density of 
the order of (10 16 GeU) 4 and negative pressure which can 
be described by an effective (and slowly varying) cosmo- 
logical parameter. Using the standard slow roll approxi- 
mation for the inflaton, the metric during the inflationary 
era can be approximated very well with a de Sitter one. 
The order of magnitude of the cosmological constant dur- 
ing the inflationary epoch is A w 1CU 11 (in Planck units). 
If other massive scalar fields are present during inflation 
such as the waterfall-field in hybrid inflationary scenarios 
[23 . [301 ] , oscillaton production is expected to be signifi- 
cant. We find that oscillatons may survive on cosmo- 
logical time-scales only if the mass of the scalar field is 
larger than 10 15 GeV/c 2 . The mass loss of an oscillaton 
during Hubble time is completely determined by the ratio 
Me/^Aj where /ia is the energy density due to A and fi c 
is that of the oscillaton in its center. For example when 
Me/^A ~ 55 the mass loss of such a "waterfall" oscillaton 
is of the order of 1% during Hubble time. This mass loss 
is exponentially decreasing for increasing values of the 
ratio /i c /ma- 

Another interesting scenario for oscillaton production 
is during reheating, assuming that the inflaton field un- 
dergoes oscillations near the minimum of its potential. 
These oscillatons store a considerable amount of energy 
and due to parametric resonances induced by quantum 
effects (for oscillons this has been shown in Ref. [24]) they 
may influence the fluctuation spectrum. 

The present value of A s=s 10 -122 (in Planck units) be- 
ing extremely small, its effect on the mass loss of generic 
oscillatons is completely negligible on a cosmological time 



scale. For example if the mass of the scalar field is of the 
order of m fa 10~ 15 eV and the amplitude of the oscillaton 
is approximatively 10 -17 , the effect of the cosmological 
constant results in a roughly 1% mass loss during the 
lifetime of the Universe. 

We also consider (spherically symmetric) boson stars 
[3ll - l34| in this setting, i.e. when a complex scalar field is 
coupled to gravity with a nonzero cosmological constant. 
Adopting our small amplitude expansion techniques de- 
veloped for oscillons and oscillatons we obtain boson star 
solutions in a systematic way. Remarkably, this method 
yields very good results even for near maximal mass bo- 
son stars. Spatially well localized boson star-type solu- 
tions with a static metric, familiar in the case A < no 
longer exist when A > 0. There exist, however, radi- 
ating boson stars with a time dependent metric loosing 
continually their mass due to emission of scalar waves. 
We obtain their mass loss using our methods developed 
for oscillons resp. oscillatons, with the result that this 
mass loss is exponentially small and it is negligible for 
sufficiently small values of A. 



II. FIELD EQUATIONS 

We consider a real or complex scalar field coupled to 
gravity. The total Lagrangian density is 



l07T 



(1) 



where the Einstein Lagrangian density in case of a cos- 
mological constant A is 



C G = 



-g(R -2A). 



(2) 



We work in a D + 1 dimensional curved space-time with 
metric g a b and assume metric signature ( — h + + ...). Un- 
less explicitly noted, we use Planck units G = c = h = 1. 
In some cases it is advantageous to use natural units, 
where c = Ti = 1 (but G ^ 1) and all quantities are mea- 
sured in various powers of electron volts. Because of the 
different notation in the potential term we consider the 
real and the complex field case separately. Variation of 
the action (fTJ) with respect to g ab yields Einstein equa- 
tions, 



Gab + A.g a h = 8nT a b 



Boson stars 



(3) 



We consider a complex scalar field $ with a self- 
interaction potential [/($*$). For a free field with mass 
to the potential is [/($*$) = m 2 $*$. The Lagrangian 
density belonging to the scalar field is 



C 



M 



[<? a6 $* a $ 6 + [/($*$)] 



(4) 



Variation of the action (TTJ) with respect to $* yields the 
wave equation 



S a& $ ;ab -$£/'($*$) = 0, 



(5) 



where the prime denotes derivative with respect to $*$. 
The stress-energy tensor in this case is 

Tab = \ {*;* 6 + ** 6 *,„ - 9ab &**%* ,d + £/($**)] } . 

(6) 
We shall assume that the self-interaction potential, 
[/"($*$), has a minimum L/ ($*<!>) = at $ = 0, and 
expand its derivative as 



£/'($*$) = m 2 



l + ^u fe ($*$) A 



fc=l 



(7) 



where to is the scalar field mass and Uk are constants. 

In order to get rid of the 87r factors in the equations 
we introduce a rescaled scalar field and potential by 



Then 



= V8tt$ 



U' 



U((t>*(t>) = 8tt£/($*$) . (8) 
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Uk( 



/.* J*\k 



fc=l 



with 



Uk 



Uk 

{8ir) k 



(9) 



(10) 



B. Oscillatons 

In this case we consider a real scalar field $ with a self- 
interaction potential U(§). For a free field with mass 
to the potential is U(<&) — to 2 $ 2 /2. The Lagrangian 
density belonging to the scalar field is 



-M 



-u i i$ a $< a + U($) 



(11) 



Variation of the action (JTJ) with respect to $ yields the 
wave equation 



g a ^;ab -[/'($) =0, 



(12) 



where the prime now denotes derivative with respect to 
<F The stress-energy tensor is 

T ab = $,a$,fc - gab U*,c*' C + U ($)) . (13) 

We shall assume that [/($) has a minimum [/(<&) = 
at $ = 0, and expand its derivative as 



[/'($) =to 2 I $ + ^M fe $* 



(14) 



k=2 



where m is the scalar field mass and Uk are constants. 

In order to get rid of the 87r factors in the equations 
we introduce a rescaled scalar field and potential by 



Then 



</> = %/8^$, U{4>) = 8ttI/($) . 



[/'((/>) = to 2 




with 



Uk 



u k 



( 87r )(fc-l)/2 



(15) 



(16) 



(17) 



C. Scaling properties 

In both the real and the complex case, if the pair 4>(x c ) 
and g a b(x c ) solves the field equations with a potential U 
and cosmological constant A, then 

4>(x c ) = ^x c ) , g ab (x c ) = g ab { 1X c ) (18) 

for any positive constant 7, is a solution with a rescaled 
potential 7 2 £7 and rescaled cosmological constant 7 2 A. 
This scaling property may be used to set the scalar field 
mass to any prescribed value, for example to make m — 1. 
If D = 1 then, by definition, the Einstein tensor is 
traceless, and from the trace of the Einstein equations 
it follows that the potential U is constant. Hence we 
assume that D > 1. 



D. Spherically symmetric D 
space-time 



1 dimensional 



We consider a spherically symmetric D + 1 dimen- 
sional space-time using isotropic coordinates x^ — 
(t, r, 9±, ..., #d_i). In this coordinate system the metric is 
diagonal and its spatial part is in conformally flat form, 
with components 



9tt = -A , g rr = B, 



n-l 



= r 2 B 



g BnBn = r 2 B J] sin 2 9 k , 
fc=i 



where A and B are functions of temporal coordinate t 
and radial coordinate r. In case of oscillatons, the pe- 
riodic change in the acceleration of the constant radius 
observers is much smaller in isotropic coordinates than in 
the more commonly used Schwarzschild coordinate sys- 
tem [III |35|. It is also more convenient to study the 
Newtonian limit of boson stars using isotropic coordi- 



nates [36] . The components of the Einstein tensor are 
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2Ai \rA?BJ , r 
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rB? 



(21) 
(22) 

(23) 



The potential independent term in the wave equation ([5]) 
and dni is 



g ab ^ ab 



j>.rr _ <j\tt _ </>,« ( B 

B A 2B D I A 



<\>,r 



2D-2 atjD-2\ 



2r 



2D-2 



AB D 



— (r 2D - 2 AB 



(24) 



E. de Sitter space-time 

If the cosmological constant A is positive, far from the 
central boson star or oscillaton the space-time should ap- 
proach the dc Sitter metric. In static Schwarzschild co- 
ordinates on a D + 1 dimensional space-time this metric 
can be written as 



ds 2 



'I - H 2 f 2 ) dt 2 



df 2 



1 - H 2 f 2 



+ f 2 dn 2 D _ 1 , (25) 



where dVt 2 D _ l is the metric on the unit sphere S D ~ l . 
The Hubble constant H is related to the cosmological 
constant A by 



H 2 = 



2A 



D(D-l) 



(26) 



In isotropic coordinates the de Sitter metric takes the 
form 



ds 2 = - 



1 



H 1 



■dt 2 



dr 2 + r 2 dn|,_ 1 



1 



H 2 r 2 



(l+^ 



, (27) 



4 / V 4 

where the two radial coordinates are related by 



1 



H 2 r 2 



(28) 



The cosmological horizon, which is at f = 1/H in the 
Schwarzschild coordinate system, is at r^ = 2/H in 
isotropic coordinates. 



III. SMALL-AMPLITUDE EXPANSION 

The small-amplitude expansion procedure has been ap- 
plied successfully to describe the core region of flat back- 
ground oscillons [l], |j, |37|, |38[ • The method was general- 
ized in [35[ to the case when the scalar field is coupled 
to Einstein gravity, and in 39] to a very similarly be- 
having scalar-dilaton system. In this section we expand 
oscillatons in the case when the cosmological constant is 
positive. 

Newtonian boson stars have already been investigated 
in the paper of Rufhni and S. Bonazzola [32j . obtaining a 
system of two coupled differential equations, which nowa- 
days are generally called Schrodinger-Newton equations 
in the literature |40l443l | . The same system was obtained 
as weak gr avity limit of general relativistic boson stars 
in [44( and [45[ . The leading order results of the small- 
amplitude expansion of boson stars, that we present in 
this section, when restricted for A = 0, also yield the 
same equations. 

The repulsive effect of the cosmological constant is 
taken into account in a similar way as in |20l-|22|. Since 
we intend to investigate localized objects, where at large 
distances the scalar field is negligible, by the generalized 
Birkhoff's theorem |46l44a |. the metric should approach 
the Schwarzschild-de Sitter metric. This also means that 
in the asymptotic de Sitter region the expansion rate 
(Hubble parameter) is constant. 



A. Choice of coordinates 

We are looking for spatially localized bounded solu- 
tions of the field equations for which cf> is small and the 
metric is close to de Sitter. The smaller the amplitude 
of a boson star or an oscillaton is, the larger its spatial 
extent becomes. Therefore, we introduce a new radial 
coordinate p by 



p — emr . 



(29) 



where e denotes the small-amplitude parameter. Moti- 
vated by the scaling property (fT5)l we have also included 
a m factor into the definition of p. 

Long-living localized solutions are expected to exist 
only if their size is considerably smaller than the distance 
to the cosmological horizon. We introduce a rescaled 
Hubble constant h by 



H = e rah . 



(30) 



and assume that even h is reasonably small. This way we 
ensure that the typical size of boson stars or oscillatons, 
which is of the order l/(em) using the r coordinates, 



remains smaller than the radius of the cosmological hori- 
zon, r h = 2/H = 2/(e 2 mh). 

We will see in Section lill El that the energy density p of 
the scalar field at the central region of small amplitude 
oscillons or boson stars is proportional to e 4 . On the 
other hand, the cosmological constant can be interpreted 
as a fluid with energy density, 



PA 



A 

8^ 



(31) 



and pressure pa = — PA- Using (j2"6")l and ([5T)]) p\ can be 
written as 



MA 



DP - 1). 4. ..-•,.■> 

16tt 



Wti 2 



(32) 



The smallness of h guaranties that p\ remains smaller 
than energy density p of the scalar field even in the small- 
amplitude limit. 

We expand (f> and the metric functions in powers of e 
as 



E 



e 2k <j) 2 k ■ 



fc=i 



A = 



B = 



1 4 



e 2 h 2 p 2 



f 



E^ 4 

fe=i 

oo 

J2 £ 2k B: 



fc=l 



(33) 

(34) 
(35) 



The first terms in A and B represent the de Sitter back- 
ground, according to (|27[). Since we intend to use asymp- 
totically de Sitter coordinates, we look for functions 4>2ki 
A2k and i?2fc that tend to zero when p is large. One 
could initially include odd powers of e into the expan- 
sions ([3"3"1> - (1531) . However, for oscillatons it can be shown 
by the method presented below, that the coefficients of 
those terms necessarily vanish when we are looking for 
configurations that remain bounded in time. It can also 
be verified directly, that the small-amplitude expansion 
of boson stars do not include odd powers of e. This is in 
contrast to flat background oscillons, where the leading 
order behavior of the amplitude is proportional to e. 

The frequency of boson stars and oscillatons also de- 
pends on their amplitude. The smaller the amplitude 
is, the closer the frequency becomes to the threshold m. 
Hence we introduce a rescaled time coordinate r by 

T = ut, (36) 

and expand the square of the e dependent factor ui as 



2 2 

lo = m 



i + Y^ e2kuj2k 



(37) 



fc=i 



It is possible to allow odd powers of e into the expan- 
sion of cl> 2 , but the coefficients of those terms turn out 



to be zero when solving the equations arising from the 
small-amplitude expansion. There is a considerable free- 
dom in choosing different parametrizations of the small- 
amplitude states, changing the actual form of the func- 
tion uj. The physical parameter is not e but the frequency 
of the periodic states that will be given by ui. Similarly 
to the asymptotically flat case in [351 ] and to the dilaton 
model in [39], it turns out that for spatial dimensions 
2 < D < 6 the parametrization of the small-amplitude 
states can be fixed by setting u — my/l — e 2 . 

The field equations we intend to solve using the r and 
p coordinates are the Einstein equations (0) , substituting 
the Einstein tensor components from (|20|) - ([23|) . together 
with the wave equation (0 or (p~2|) . We note that these 
equations are not independent. The (r, p) component of 
the Einstein equations is a constraint, and the wave equa- 
tion is a consequence of the contracted Bianchi identity. 



B. Boson stars 

In case of boson stars the metric is static, hence A}. 
and Bfc are time independent. The complex scalar field 
4> oscillates with frequency u>, 



(j) = ipe % 



(38) 



where ip is real and depends only on the radial coordinate 
p. The small-amplitude expansion coefficients of <f> satisfy 



4>k = ipke* 



4>ke v 



(39) 



where if)]- are real functions of p. Substituting the expan- 
sion (|33 |) - (|35l) into the field equations, if D ^ 2, then to 
leading e 4 order the (p, p) component gives 



B 2 = 



-4, 



2- D 



(40) 



Since for D — 2 there is no solution representing a local- 
ized object we assume D > 2. Then the (r, r) component 
and the wave equation yield to leading order a coupled 
system of differential equations for A 2 and tp 2 > 



dp 2 



D - 1 dA 2 
dp 



D 



f> 



D- 1 



1%, 



d 2 i> 2 D - 1 # 2 



dp 2 



P 



^ = "02(^2 - W 2 - h 2 p 2 ) . 

dp 



(41) 
(42) 



Introducing the functions s and 5* by 

s = lu 2 - A 2 , S = i/j 2 \2 



D-2 



D- 1 
equations (14"TT) and (l4"2"j) can be written into the form 

d 2 S D-l dS 
"dp 2 "*" 



(43) 



p dp 
d 2 s D 

dp^ + ~ 



+ (s + h z p z )S = 



(44) 



-J-1± + S 2 = 0. (45) 

p dp 



Since the same equations will be obtained for small- 
amplitude oscillatons in the next subsection, we postpone 
the discussion of this system to Subsection MIDI 

Proceeding to higher orders, the e 6 components give a 
system of lengthy differential equations linear in A±, -B4, 
ipi and their derivatives. These equations have nonlinear 
source terms containing A 2 , B 2 and ip 2 . They will also 
involve the first coefficient u\ from the series expansion 
of the potential U. To leading order, the structure of 
small-amplitude boson stars only depends on the mass 
parameter m. 



C. Oscillatons 

In this case the scalar field is real, and the metric com- 
ponents are also oscillating. We are looking for solutions 
that remain bounded as time passes. It turns out, that 
these configurations are necessarily periodically oscillat- 
ing in time. 

From the e 2 components of the field equations follows 
that B 2 only depends on p and that 



<t>2 = P2 cos(r + 5) , 



(46) 



where two new functions, p 2 and S are introduced, de- 
pending only on p. From the e 4 part of the (r, p) com- 
ponent of the Einstein equations it follows that B4 can 
remain bounded as time passes only if 5 is a constant. 
Then by a shift in the time coordinate we set 



5 = 0. 



(47) 



This shows that the scalar field oscillates simultaneously, 
with the same phase at all radii. 

From the e 4 component of the difference of the (p, p) 
and (#1 , 9i ) Einstein equations follows that A 2 also de- 
pends only on the radial coordinate p, showing that the 
metric is static to order e 2 . This is the main advantage 
of using isotropic coordinates over Schwarzschild coordi- 
nates. But even in this system the coefficients A4 and 
B4 will already contain cos(2r) terms. From the e 4 com- 
ponent of the (r, p) Einstein equations follows that 



B d 



Pi 



4(D-1 



■ cos(2t) , 



(48) 



where 64 is a function of p. The (p, p) component shows 
that if D ^ 2 then 



So 



A-2 



2- D 



(49) 



If D = 2 then A 2 = 0, and there are no nontrivial local- 
ized regular solutions for B 2 and p 2 , so we assume D > 2 
from now. The (r, r) component yields 



d 2 A 2 
dp 2 



D 



1 dA 2 
dp 



D 



D-l 



Pi- 



rn 



The wave equation provides an equation that determines 
the time dependence of (j>^. In order to keep c/>4 bounded, 
the resonance terms proportional to cost must vanish, 
yielding 



d 2 p 2 , D - 1 dpi 

-—=p 2 (A 2 



dp 2 



P dp 



LU 2 



h 2 P 2 ) 



(51) 



Equations ([50]) and ([ST]) do not depend on the coef- 
ficients Uk of the potential U(<fi). This means that the 
leading order small-amplitude behavior of oscillatons is 
always the same as for the Klein-Gordon case. 

Introducing the functions s and S by 



co 2 



S =P2 



D 



D-l ' 



(52) 



equations (|50|) and (J51I) can be written into the form 
(|44|) and (j45|) already obtained at the small-amplitude 
expansion of boson stars. 



D. Schrodinger-Newton equations 

The equations describing both small-amplitude bo- 
son stars and oscillatons on an expanding background 
has been written into a form, which for h = reduces 
to the time-independent Schrodinger-Newton (SN) equa- 
tions S3-S1 



d 2 S D-ldS . u2 3 , n 
-r-2 + i- + (s + h 2 p 2 )S = , 



dp 2 p dp 

d 2 s D-l ds 
dp 2 p dp 



S 2 = . 



(53) 
(54) 



If S(p) and s(p) are solutions of (|S"3"|) and (|5"4")) . then 
the transformed functions 

S(p) = X 2 S(\p) , s(p) = A 2 s(Ap) (55) 

solve the equations 

d 2 S D-l dS 



dp 2 



+ {s + h 2 p 2 )S = 0, (56) 

^ + ^£ + # = 0, (57) 



p dp 
d : 



which are obtained from the SN equations by replacing 
the constant h by h = X 2 h. 

If the scalar field, and consequently S, tends to zero 
for p — > oo then (|54[) gives the following approximation 
for s. 



s = s + sip 



„2--D 



(58) 



Since we are looking for solutions for which A 2 = lu 2 — 8 
also approaches zero for large p, necessarily us 2 — sq. In 
order to make the solution of the SN equations (|53l) - (|5"4")) 
unique we demand 



^2 = SO 



(59) 



Setting Wfe = for k > 3 this also fixes the connection 
between the frequency uj and the small-amplitude param- 
eter e, giving 



= my 1 



cj = my 1 — e 
Substituting into (|53[) we get 
d 2 S D-l dS ,,22 



(60) 



dp 2 p dp 



+ (h 2 p 2 -l + s lP 2 - u )S = 0, (61) 



determining the large p behavior of S. 



E. Energy density 

The energy density of the scalar field is p = T a i,u a u b , 
where the unit timelike vector u a has the components 
(l/vJ4,0, ...,0). In terms of the rescaled scalar field </>, 
the energy density of oscillatons can be written as 



li = 



1 

167T 



-AtPUt) 2 ^ 



and for boson stars 



P 



1 


[1#*# l#*^ rr,j.*^] 





— — - — f H 7— -r 1 + U(6*d>) 


16tt 


[A dt dt B dr dr yv v '\ 



(62) 



(63) 



In case of boson stars, using (|33j) . (1391) and (|43j) , to lead- 
ing order in e we get 



1 



,13 — 1 



/' 



16tt D - 2 



-m 2 S 2 



(64) 



For oscillatons we obtain the same expression for p using 
(1331) , (|4"6")l and (|52l) . Measuring the scalar field mass m in 
eV^/c 2 units, the energy density at the symmetry center 
for D = 3 can be written as 



p c 



mc 

~eV 



1.436 x 10 



39 ^9 

m ! 



(65) 



(We use Roman m for meters in order to distinguish from 
the scalar field mass to.) Here we have used the numer- 
ically obtained central value of S, which is S c — 1.0215 
for D = 3. 

The effective energy density pa corresponding to the 
present value of the cosmological constant is given in or- 
dinary units by (|D3I) in Appendix [D] Comparing to (|65l) . 
one can see that the energy density of the scalar field is 
much higher than p\ unless e or to is extremely small. 
Using (j32|) . the ratio of p\ to the central density p c can 
be written as 



Pa 



D{D-2) 

S 2 



h 2 



(66) 



providing an important physical interpretation for the 
rescaled Hubble constant h. In Table U we give the value 
of S c and the coefficient of h 2 in pa/Pc for the relevant 
dimensions. We will see in the following sections that h 
turns out to be the essential parameter determining the 
energy loss rate of boson stars and oscillatons. 





L> = 3 


D = 4 


D = 5 


s c 

D{D-2)/S 2 c 


1.0215 
2.8751 


3.5421 
0.63761 


14.020 
0.076313 



TABLE I: The numerically calculated values of the function 
S at the center and the coefficient of h in /Ua//Uc for D = 3, 
4 and 5 spatial dimensions. 



F. Outer core region 

If 2 < D < 6, and the cosmological constant is zero, 
i.e. /i = 0, then assuming that so = — 1 the SN equations 
have a unique localized nodeless solution. Solutions with 
nodes have higher energy and are unstable. For large p 
the function S satisfies the equation 



d 2 S 
dp 2 



D - 1 dS . 

-7- + (-l 

p dp 



s\p 



„2-D 



)S = 0, 



which has exponentially decaying solutions 

T 



S = S t 



S = S t 



,l-«i/2 



p (D-l)/2 



1 + 



P 

oil 

\P 



D 



D>3. 



(67) 

(68) 
(69) 



The numerically determined values of the constants Si 
and St for the case h = are given in Table |TIJ 





D = 3 


D = 4 


D = 5 


Si 

Si 


3.50533 
3.49513 


7.69489 
88.2419 


10.4038 
23.3875 



where a and /3 are constants, describing the amplitude 
and phase. We can interpret this standing wave as the 
superposition of an outgoing wave carrying energy out 
from the central object and an artificial ingoing wave 
component that is added to keep the solution exactly 
periodic. We will relate the amplitude of the tail to the 
energy loss rate of dynamical boson stars or oscillatons. 
When solving the SN equations (JSHJ) and (|M|) by a nu- 
merical method, the amplitude a depends on the central 
values of s and S given at p = 0. We intend to determine 
the minimal tail amplitude, a = a m - ln , assuming that s 
tends to So = — 1 for p — > oo. We are interested in solu- 
tions which has no nodes in the central region where the 
influence of h is negligible. If h > then a m i n is nonzero. 
In order to show the influence of h on the form of the 
function S, on Fig. Q] we plot the radial dependence of S 
for h = 0.1016 and for h = 0. 
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TABLE II: The numerical values of the constants si and St 
for h = in 3, 4 and 5 spatial dimensions. 

Even if h is nonzero, we assume that it is small enough 
such that there is a region of p which is well outside the 
core region, but where the influence of the cosmological 
constant is still negligible. In this region S satisfies (I571) . 
and consequently (|68[) and (|69p are good approximations 
together with the values of the constants belonging to 
the h = case given in Table [TTJ 



G. Oscillating tail region 

If h is nonzero, at very large distances p 3> l//i, and 
the behavior of S is determined by the equation 



d 2 S D-l dS 



dp 2 p dp 



h'p'S = . 



(70) 



The function S has an oscillating standing wave tail in 
this domain, 



S = -^ cos (-£- + f3 



(71) 



FIG. 1: The radial dependence of the function S for a large 
value of h is compared to the h = flat background case. The 
function S for h = decays exponentially, while for h > 
curve starts to oscillate. For h — 0.1016 the first change in 
signature is at p t — 8.04. 

Assuming that h is small, we calculate the amplitude 
a by WKB analysis. Our aim is to continue the function 
S given by (|68|) and (|69|) into the very large distances 
where the cosmological term h 2 p 2 dominates. We will 
express the amplitude a in terms of the constant St ■ To 
the function which is exponentially decaying when going 
outwards at the outer core region belongs an oscillating 
tail with a certain amplitude and phase. Adding the 
mode which exponentially grows when going outwards 
does not change the function noticeably in the outer core 
region but adds an oscillating tail with a phase shift 7r/2. 
The resulting tail amplitude is minimal if this second 
mode is not present. Hence the WKB method will give 
the minimum amplitude a = a min . In Sec. IIV Fl we will 
compare the WKB amplitude to results obtained by di- 
rect numerical integration of the SN equations. Besides 
indicating the correctness of the WKB results, this will 
also give information on how large h can be chosen such 
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that the WKB result can be still considered a valid ap- 
proximation. 



IV. WKB FORMALISM 



A. Schrodinger form 

Introducing a new radial coordinate 
V = hp 



(72) 



we transform (|61[) into a form corresponding to a 
Schrodinger equation with h playing the role of the 
Planck's constant Ti, 



, , . d 2 S D - 1 dS ■ , , 
\dy 2 y dy ' 



,D-2 



y 



D-2' 



Sl 5 = 0. 



Introducing 

we can transform this into a one-dimensional form, 
h 2 ^-+p 2 Z = 0, 



(73) 
(74) 

(75) 



where 



2 2 1 

P =y - 1 



yD-2 4^2 



(D-l)(D-3). (76) 



The outer domain of the core region is represented now 
by < y <C 1, where 



Z = S t h l - Sl / 2 y Sl / 2 e-' 


l/h 


i + o( 


7/ / 


, D = 3, (77) 


z = s t h^ D -^ 2 e-y/ h 


14 


*G): 


) 


D > 3 . (78) 



B. WKB approximation 

Introducing a new function / by 
Z = e lf/h 
and substituting into (fT5"j) we get 

d 2 / /d/ 



dy 2 Vdy/ 
We write f as a power series in h, 

oo 

fc=0 



(79) 



(80) 



(81) 



Differently from the standard WKB approach, now p 2 
also has to be expanded in h, 



v 2 = Y.v\h\ 



fe=0 



where the nonvanishing coefficients are 

po = y 2 - 1 , 

2 _ (D-l)(£>-3) 

P2 — ~~ 



V 



PD-2 ~ n_ 2 



Substituting into (I5D1 . to leading order we get 



dfo\ 2 d fo , 



The order /i components yield 

i fJo + 2 = 2 dhdh 
dy 2 dy dy 



(82) 

(83) 
(84) 
(85) 

(86) 
(87) 



C. WKB for D > 3 

If D > 3 then p 2 = 0, 

and we get the standard WKB result, 

Z = -^= exp (±~ f y pody) , < y < 1 , (89) 
VPo V h J\ J 



Ipo 

Z = ^expi+-J lh) ,hi). ,, - I. (<)()) 



The general solution for < y < 1 is a sum of two terms 
with proportionality constants A+ and A_ , and similarly 
for the y > 1 case with constants B + and _B_. For the 
< y < 1 case |po| = \/l — y 2 - Substituting (|83|) . it is 
possible to perform the integral in the exponentials. For 
0<y<l, 



hh 



Pody 



hh 



\A - V 2 dy 



'2h( V ^y 2 



arcsnn/ - -j (91) 



while for y > 1, 

i rv { rv 

r / Po^y = r / v y 2 - i rf y 



2h 



i\fy 2 - i - in ( y + Vy 2 - i 



(92) 



In the region < y -C 1, substituting the expansion of 
nj, to leading order ([B"9")) takes the form 



Z = A± exp 



T 



1 / 7T 



(93) 



Comparing with ([78)) , follows that for D > 3 

A+ = S^- 1 )/ 2 exp (-^) , A_ = . (94) 



D. WKB for D = 3 

Because of the term proportional to h in p 2 , the I? = 3 
case has to be treated separately. Then ([8"7| takes the 
form 



i d n co. m a. Sl d /i 



Using 



(95) 



T- ■ ( 96 ) 



i d 


In f+t/i 


/2-1) 


1 1 


dy 


2dyV \ va I 


J 2tyVl - 2/ 2 


Integrating, for < y < 1, 


/ 1 = iln(Vj/ 2 - 







2i 


lny — 


n (l + Vi - y 2 ) 


+ Ci, 


and for j/ > 1, 






'-•) 


=F — arct; 
2 


/ ' \ 

in 1 1 -(- 


C2- 


w^-v 


Substituting into ([75]). for < y < 1, 


A 
£= r- 

V^ 


/ \ ±S1 ' 2 / 


^V?, 


po| Vl + Vl-2/V ^ ; 


< 
and for y > 1, 


* = /— ex P < T 


f i 


•)-; 


1 


Vx/s 2 -l 


J 



(97) 



(98) 



(99) 



X 6XP \ /i / P ° dy 
For < y <C 1 the expansion of (|9"9"|) gives 

£ = A± (JJ exp 

Comparing with (|77|l . follows that for D = 3 
A+ = 2Sl /2 5 ^i- Sl /2 exp / ^L\ A_=0. (102) 



T 7^~4 



(101) 



E. Connection formulae 

The next step is to relate the amplitudes A± in the y < 
1 region to the amplitudes B± for y > 1. The expressions 
(J591 and flU)) for D > 3 are just the standard WKB 
solutions, so one can apply the appropriate connection 
formulae to express B± in terms of A±. In our case we 
need the one where the field is exponentially increases 
when goin g aw ay from the turning point (see for example 
(34.18) of ggj]) 



:COS 



Pody 



1 

1 

-, exp . 



bolrfy 



from y > 1 to y < 1. This is equivalent to 

A+ ( in 

A_=0, J B± = ^texpf± T 



(103) 



(104) 



The expressions (|99|) and (|100p for I? = 3 differ from the 
standard WKB solutions by factors involving si. How- 
ever, since both of these factors take the value 1 at the 
y = I turning point, the formulae (|104[) relating the am- 
plitudes at the two sides hold in this case too. A detailed 
derivation of (|104j) is given in Appendix [A] 



F. Schrodinger-Newton tail amplitude 

For y ^> 1 the expression (|92[) can be approximated by 

(105) 



i f y , ^ 2 



2/i 



In the D > 3 case to leading order (IM1) yields 



Z = exp ± — 



Substituting B± from (|104[) gives 

v A + (V* .* 

Z = —P cos \TT + T 

V^ V2/l 4 

For D — 3, to leading order (|100[) yields 



£>> 3 






exp 



±; 'lf 



7TSl 



V^O) an( j substituting B± from (|104[) gives 



Z = — cos 



§; + !<« + » 



D 



(106) 



(107) 



(108) 



(109) 



Using ([74]) we can express the original S variable in the 
SN equations using the coordinate p = y/h, obtaining the 
same asymptotic formula as in (|71[) . 



S = -^ cos (-£- + p 



(110) 
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where the amplitude is 



h D/2 



(111) 



for any D 1 and the phase is 



P={ 



IT 

4 

7T 
I 4 



(si + 1) if£> = 3 
if L>>3 



(112) 



Substituting the value of A + from (|94| and (|102|) gives 
our final WKB result for the minimal amplitude 
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FIG. 2: The curves represent the WKB tail amplitude, while 
(113) the points the numerically calculated values, as a function of 
h, for D = 3, 4, 5 spatial dimensions. 



It is instructive to compare the above analytical result 
to the minimal amplitude obtained by the numerical so- 
lution of the SN equations (l53l) - (l54l) . For a chosen h one 
has to minimize the oscillating tail for those solutions for 
which the field s tends to sq — — 1 asymptotically, while 
S has no nodes in the core region. The solutions depend 
very strongly on the central values s c and S c . The results 
for various h for D — 3 spatial dimensions are listed in 
Table IIIII The amplitudes ct m in for D = 3, 4, 5 dimen- 



h 


8c/8t° 


Sc/SZ=° 


Pt 


„ nura 


"mm 1 


Q WKB L 


0.036111 


0.96295 


0.97043 


28.65 


7.7377 x 10" 6 


0.033 


0.046769 


0.93539 


0.94830 


22.04 


6.2645 x 10" 4 


0.054 


0.064593 


0.86123 


0.88802 


15.73 


0.033101 


0.120 


0.071112 


0.81742 


0.85171 


14.13 


0.085425 


0.177 


0.076502 


0.76580 


0.80827 


12.90 


0.16998 


0.268 


0.085350 


0.63180 


0.69420 


10.62 


0.49007 


0.614 


0.088000 


0.58803 


0.65720 


9.94 


0.64553 


0.727 


0.091400 


0.53709 


0.61443 


9.24 


0.85085 


0.778 


0.094800 


0.49518 


0.57976 


8.72 


1.02507 


0.709 


0.101600 


0.43960 


0.53576 


8.04 


1.25256 


0.402 



TABLE III: For several choices of h the numerically calculated 
minimal tail amplitude aJS™ and its relative difference from 
the WKB value (|113l) is listed assuming D — 3. We also 
give the ratio of central values of the functions s and S to 
the flat background values s h c =0 = 0.93832284 and S^ =0 = 
1.02149304. In order to show where the tail region starts, we 
also list the radius pt, where S first crosses the value zero. 

sions are plotted on Fig. and the relative difference 
from the analytically obtained WKB value is displayed 
on Fig. [3J From this figure we read off that the WKB 
analysis gives less than 10% error, if h < 0.06 in case of 
D = 3, and if h < 0.1 when D = 4 or 5. 
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FIG. 3: Relative difference of the numerically calculated min- 
imal tail amplitude from the analytically calculated WKB 
values. 



G. Scalar field tail amplitude 



For oscillatons, using (|33J}, (gBJ), f>2) and flTTO")) . the 
real scalar in the tail region behaves like 



= e 2 (j)2 = e 2 p 2 cost = e 2 S 



D-l 

D-2 



■ COST 



j) A I lit mr 

^cosM)cos( !- )) . 11141 



where the amplitude is 



e 2 



m 2 



% V D-2 



(115) 
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Applying the identity 
1 



cos a cos b = - [cos(a — b) + cos(a + b)] 



(116) 



it is apparent that (|114[) is a sum of an ingoing and out- 
going wave, both with amplitude (J>a/2. One might be 
tempted to add an ingoing wave with an opposite am- 
plitude to obtain a purely outgoing solution, and to con- 
clude that the remaining outgoing component has the 
amplitude 4>aJ2. However, an ingoing wave contains 
a cos (he 2 m 2 r 2 /2 + 0) component, which would change 
the amplitude of the WKB mode which is increasing 
when going away from the turning point. What one 
can do instead is adding a standing wave proportional 
to sin (he 2 m 2 r 2 /2 + /3) , corresponding to the suppressed 
WKB mode. If the time dependence is also phase shifted, 
then we obtain the minimal amplitude outgoing wave, 



he 2 m 2 r 2 



+ j3-ut 



In case of boson stars, using (l33l) . (j39|) . (J43J) and (jllQI) . 
the complex scalar in the tail region behaves like 



= £ 2 </> 2 = e 2 4' 2 e' T = £■ 



S D-l 



V2 



A e^'cos 



V2V D-2 

? 2 

m r 



■0 



(118) 



Now one can obtain the minimal amplitude outgoing 
wave by adding a term 



-2^2^,2 



,_|^ 31 „(^ : 



obtaining 



V2 



— exp 

f 2 



( he 2 m 2 r 2 



+ 0-ut 



(119) 



(120) 



H. Mass loss rate 



Even if p 3> l/h in the oscillating tail region, according 
to (|M|) and (|55|) , for small e the metric functions A and B 
are still very close to 1 . The mass loss rate of the system 
can be calculated by applying the expression (|Blip for 
the time derivative of the mass function from Appendix 

m 



dm 
~~dt 



2tt 2 r 



M-^D-X 



r(#) 



-T t , 



(121) 



Substituting the stress energy tensor from 

dm ir^~r D ~ 1 
~dT 



or £13 



r(#; 



($* * p + $* r $ t ) 



sr(f) 







(122) 



which expression is valid for both real and complex fields. 
We have to evaluate this expression at large r in the tail 
region. Substituting the outgoing wave form (|117|) and 
averaging for one oscillation period, or substituting (I120p . 
for both the oscillaton and boson star case we obtain 



dm 
~dT 



TT 2 



~ 1 he 2 m 3 (j) 2 A 



sr(f) 



(123) 



Here, since we are interested in the leading e order result, 
we substituted lo = m, and for large r we neglected the 
term arising when taking the derivative of the r~ D ' 2 fac- 
tor. For large radius m agrees with the total mass M of 
the oscillon or boson star. Using (|1 13|) and (|1 15|) for the 
mass loss rate and substituting h — H/(me 2 ), for D > 3 
we obtain 



dM 

~~dT 



HD 



l)S 2 e e - D 



8m^- 3 ( J D-2)r(f) 



exp 



Time 
2H 



(117) if £) = 3 there is an extra term involving s\, 



dM 
~dT 



Sfe 3 



/2me 2 



exp 



Time 
2H 



(124) 



(125) 



For small-amplitude configurations the parameter e 
can be expressed by the mass M using (|B14[) . Because 
of the complexity of the resulting expression we only give 
the result for D — 3, when M = esi/(2m) 



dM 
~dT 



-45? 



m 3 M 3 /8m 3 M 2 



Hs 2 



exp 



2-Km z M 2 
Hsf~ 



(126) 

For specific m and H this can be integrated numerically 
in order to obtain the change of the mass as a function 
of time. 

Using the leading order expression (|B14[) for the total 
mass and dividing by M we can obtain the the relative 
mass loss rate. The Hubble time Th = 1/H describes the 
time-scale of the expansion of the universe. The value of 
Tjj obtained from the present value of the cosmological 
constant is given by (ID7[) in Appendix [Dj Dividing by 
both M and H we obtain the relative mass loss extrapo- 
lated for a period corresponding to the Hubble time. For 
D > 3, 



T H dM 

~M~d7 

and for D = 3, 

T H dM 
~M~dT 



S 2 



(D-2)sih 



S 2 t 



exp 



2h 



s\h 



V exp 



V 2h) ' 



(127) 



(128) 



where the values of the constants s\, St are given in Ta- 
ble Ql] This result is independent on the units used to 
measure the time and mass, and it depends on H, m 
and e only through the rescaled cosmological constant 
h = H/(me 2 ). According to (|66| the square of the pa- 
rameter h is proportional to the ratio of the central den- 
sity /i c to the effective energy density corresponding to 
the cosmological constant /ia. On fig. [4] and on Table HVl 
we give j*j- 4M as a function of the density ratio — . 
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large e the expression (|129|) can be applied. Apart from 
the amplitude parameter e the mass loss (|125l) depends 
only on the ratio H/m. It can be checked numerically 
that for a concrete choice of H/m there is only one e 
value where the two type of mass loss rates are equal, 
which we denote by e e . For low H/m values the function 
e e can be expanded as 



C3 



with 



si 



H 



R (H 
3cJ \m 



1/3 



(l + ei+ea + ...) , 



1/3 



£2 



C4R 



H\ 2/3 
m J 



(131) 



(132) 



FIG. 4: The relative mass loss during a Hubble time period, 
-»f-^ff) as a function of — for D — 3, 4 and 5 spatial dimen- 



T H dM 
M dt 


L> = 3 


D = 4 


D = 5 


h 


Me 

MA 


h 


Me 

MA 


h 


Mc 

MA 


1 

0.1 
10~ 2 
10" 3 

10- 4 

10" 5 


0.11807 
0.09460 
0.07985 
0.06951 
0.06176 
0.05569 


24.95 
38.87 
54.55 
71.99 
91.19 
112.1 


0.20053 
0.15073 
0.12139 
0.10188 
0.08791 
0.07739 


39.00 
69.03 
106.4 
151.1 
202.9 
261.9 


0.42136 
0.23793 
0.16998 
0.13327 
0.10999 
0.09382 


73.81 
231.5 
453.5 
737.8 
1083 
1489 



TABLE IV: The rescaled cosmological constant h and the 
density ratio — belonj 
a Hubble time period. 



density ratio — belonging to various relative mass loss during 



V. MASS LOSS RATE OF OSCILLATONS IN 
THE A = CASE 

If the cosmological constant is zero and the space-time 
is asymptotically flat then the metric of boson stars is 
static and their mass is independent of time. In contrast, 
oscillatons are very slowly losing mass by scalar field ra- 
diation even in the A = case [16| . This mass loss is 
nonperturbatively small in the amplitude parameter e. 
According to [331, tor D = 3 this mass loss rate depends 
only on the amplitude parameter, and it is given by 



dM 

~dT 



Cl 

— exp 



(-7) 



e* v e 
where the numerical values of the constants are 

ci=30.0, c 2 = 22.4993. 



(129) 



(130) 



For small e this is exponentially small, while the mass 
loss induced by the cosmological constant, given by (1125[) . 
tends to zero only polynomially when e — ¥ 0. This implies 
that if A =^ then for small e the mass loss (I125P originat- 
ing from the cosmological constant dominates, while for 



■K V Cl 



,'c 4 2si\ ( H 
C4 m c 3 + [ — I In 



where the newly introduced two constants are 



2c 2 



1/3 



c 3 = — . c 4 = - (2si + 5) 

IT J IT 



m 
(133) 



(134) 



On Fig. [S]we plot e e as a function of H/m. To illustrate 
how well the above analytic approximation works we plot 
for relatively large values of H/m, although the solutions 
for e > e max w 0.5 are unstable. The quantity H/m is 




FIG. 5: The values of e e , for which the oscillaton mass loss 
rates (|125l) and (|129p are equal are plotted as a function of the 
ratio of the Hubble constant H and the scalar field mass m 
(solid line). Above the line the asymptotically flat result (|129p 
is dominant, while below (I125|l induced by the cosmological 
constant gives a larger contribution. The first three orders of 
the analytical approximation given by (|131[) are also shown. 



originally assumed to be in Planck units, but it has the 
same simple form in natural units if the Hubble constant 
is expressed in electron volts. 
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VI. POSSIBLE EFFECTS IN THE CURRENT 
UNIVERSE 

From the definition (I3U1) of the rescaled cosmologi- 
cal constant we can express the square of the small- 
amplitude parameter e. The amplitude of the scalar field 
is proportional to e 2 to leading order. Measuring the cos- 
mological constant in 1/s units and the scalar field mass 
in eV/c 2 units, (|3T))) can be written as 



, 6.58 x 1CT 16 TT eV 

£ = h HS 2 



:t3 ■ ( 135 ) 



According to (|127|) and (|128[) , the part of the mass lost 
during a Hubble time period only depends on the value of 
the rescaled Hubble constant h. For example, we can look 
for configurations where this mass loss is one percent. 
According to Table HVl for the physically interesting 3 + 1 
dimensional case this corresponds to the value of /lo.oi = 
0.07985. Using the present value of the Hubble constant, 
which is given by (ID6|i . 



fc 0.01 



w eV 
1.60 x 10- 32 — Ti- 
mer 



(136) 



Oscillatons or boson stars are generally stable if the 
amplitude parameter e, which for larger am plitudes can 
be defined from the frequency w by e = yl — us 2 /m 2 , 
is smaller than e max w 0.5. From (|136|) it is apparent 
that significant mass loss for the larger amplitude con- 
figurations with e ss £ m ax can happen only for extremely 
small scalar field masses. Very small e values are not 
excluded, but since for D = 3, according to (|B16|) . the 
total mass of the oscillaton or boson star is proportional 
to e, these configurations tend to have extremely small 
masses. Similarly, since by (|C6[) the radius is inversely 
proportional to e, small e values belong to spatially ex- 
tended configurations. It is obvious that the influence of 
the cosmological constant is larger for large radius ob- 
jects. In Table IV1 we list the values of e, the total mass 
and the radius for various choices of the scalar field mass 
m, assuming that the mass-loss extrapolated for a Hub- 
ble time period is 1%. Keeping the same scalar mass, for 
larger e, and consequently, for larger total masses and 
smaller radiuses, the mass loss rate is smaller than 1%. 
For all configurations listed in Table W\ the A = oscil- 
laton mass loss rate given by (|129l) is negligible, many 
orders of magnitude smaller than 1%. 

To show some concrete examples, we examine three 
specific choices for the scalar field mass. Firstly, for an 
axion with m = 10 _5 eV^ the mass of the oscillaton for 
which ^S-^M- — 0.01 is about that of a very small aster- 
oid, and its radius is about 15 astronomical units. Sec- 
ondly, if m — 9.55 x 10 _18 eV/c 2 , then the total mass is 
about a solar mass, M = 2.00 x I0 30 kg, which is dis- 
tributed in a region of radius of 238 light-years. Thirdly, 
for to = IO -25 the corresponding mass is about the mass 
of the Milky Way (including dark matter), and the ra- 
dius is 2 x 10 6 light-years, which is about four times the 



mc /eV 


£0.01 


M/kg 


rgs/m 


10 10 


1.27 x 


io- 21 


5.90 x 10 -11 


6.95 x 10 4 


10 5 


4.01 x 


io- 19 


1.87 x IO" 3 


2.20 x 10 7 


1 


1.27 x 


io- 16 


5.90 x 10 4 


6.95 x IO 9 


io- 5 


4.01 x 


io- 14 


1.87 x 10 12 


2.20 x 10 12 


10 -io 


1.27 x 


io- 11 


5.90 x IO 19 


6.95 x 10 14 


io- 15 


4.01 x 


10 -9 


1.87 x IO 27 


2.20 x 10 17 


1Q -20 


1.27 x 


io- 6 


5.90 x 10 34 


6.95 x IO 19 


10 -26 


4.01 x 


10" 4 


1.87 x 10 42 


2.20 x 10 22 


IO" 30 


1.27 x 


io- 1 


5.90 x 10 49 


6.95 x 10 24 



TABLE V: Oscillaton or boson star configurations for which 
1% of the mass is lost during Hubble time, assuming the 
present value of the cosmological constant. For each choice of 
the scalar field mass m the amplitude parameter e, the total 
mass M and the radius containing 95% of the mass is given. 



radius of the stellar disk. According to Table IIVI and 
(|D3[) . for all the above configurations the central density 
is fi e = 54.55/ja = 3.7 x 10- 22 kg/m 3 . 

The maximal amplitude boson stars or oscillatons, 
with e « e ma x generally have much bigger masses and 
smaller sizes than the states with 1% mass loss discussed 
in the previous paragraphs. For example, the e = 0.5 os- 
cillaton for the axion case has a mass M = 2.3 x I0 25 kg 
(about four Earth masses) and radius 7-95 = 18cm. Only 
considering the mass loss (|125[) induced by the cosmo- 
logical constant, the maximal amplitude boson stars or 
oscillatons have extremely long lifetime. For the ax- 
ion case, using (|135p . we get h — 5.1 x IO -28 , and 
the relative mass loss during a Hubble time period is 
of the order 10~ 10 , certainly negligible. For boson stars 
formed by complex fields this is the only mechanism to 
lose mass, consequently, close to maximal mass boson 
stars have practically constant mass. However oscillatons 
lose mass even in the A = case by (|129p which gives 

m"1F" = - 3 - 1 x 10l ° for £ ~ °- 5 ' lowing that this state 
is unstable on the cosmological time-scale. However, the 
mass loss decreases exponentially with s. An initially 
maximal mass oscillaton created in the early universe 
loses mass relatively quickly, and its amplitude param- 
eter decreases to about e = 0.31 during a Hubble time 
period, where according to (|129|) the relative mass loss 
rate is ^-M-^M- = —0.14. For more details and different 

M at 

scalar field masses sec Tabic VIII of 35]. 



VII. INFLATIONARY ERA 

The energy density of the inflaton field during the in- 
flationary epoch can be estimated as /j,\ = (10 16 GeF) 4 
(see e.g. [26jh This expression is valid in natural units, 
where c = h = 1 (but G^l) and the energy is measured 
in electron volts. In ordinary units this corresponds to 
the mass density /z\ = 2.3 x 10 84 fcg/m 3 . Assuming a de 
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Sitter geometry, the Hubble constant is 



//= l /^^=3.6xl0 37 i=2.4xl0 13 G e l/, (137) 
3 s 

in ordinary and natural units. By (|26[) this belongs to a 
cosmological constant 



A 



3H 2 



4.3 x 10 



5 s 



1 



1.1 x 10 



-li 



(138) 



c m' 

in ordinary and Planck units, respectively. 

If in addition to the inflaton field there is another scalar 
field x on this de Sitter background, and it is massive, 
then it is likely to form localized oscillaton configurations 
(or boson stars in case of a complex field). This second 
field may be, for example, the waterfall field in a hybrid 
inflation theory [29J, |30j. If the mass parameter of the 
field x is to, then we can apply (I135[) to relate the small- 
amplitude parameter e to the rescaled Hubble constant h. 
The parameter h is the essential parameter determining 
the mass loss rate of the object. For example, if the mass 
loss extrapolated for a Hubble time period is assumed to 
be 1%, then h = ho.oi = 0.07985, and for the inflationary 
era we obtain 

3.0 xlO 14 ^. (139) 



£ 0.01 



Since oscillatons and boson stars are stable only if the 
amplitude parameter satisfies e < £ max s=s 0.5, states 
with 1% mass loss can only exist for to > 1.2 x 
!0 15 GeV/c 2 . For objects formed from fields with to < 
1.2 x !0 15 GeV/c 2 , necessarily e < e max < e .oi- Then 
from (|135p it follows that h > /lo.oij which implies that 
the mass loss during the Hubble time is always larger 
than 1%. Oscillatons and boson stars formed from scalar 
fields with mass to < 10 15 GeV/c 2 are very short living 
during the inflationary era. Another way to see why this 
lower limit is necessary is to calculate the central den- 
sity n c by (|65[) . and observe that for small scalar field 
mass [i c < ^a, and hence long living oscillatons are not 
expected to exist. 

For scalar field masses to > 1.2 x !0 15 GeV/c 2 we list 
the properties of some oscillatons with 1% mass loss in 
Table IVII All these states have the central density /i c — 
54.55/m = (2.7 x 10 16 GeV) 4 = 1.3 x 10 86 fcg/m 3 , and the 
A = mass loss is much smaller than 1% for them. 

The mass loss rate induced by the cosmological con- 
stant decreases very quickly when considering higher am- 
plitude configurations. The behavior of the A = os- 
cillaton mass loss is just the opposite. The two kind 
of relative mass loss extrapolated for a Hubble time 
period of some maximal amplitude states belonging to 
£ = £max ~ 0.5 are given in Table IVTll 



VIII. REHEATING 

Oscillatons are likely to form after the end of inflation, 
in the reheating era, when the inflaton field is oscillat- 
ing around the potential minimum. These oscillatons are 



mc 
GeV 


£0.01 


M 
kg 


125. 
m 


M ,.™„ 


1.2 x 10 15 
10 16 
10 17 
10 18 
10 19 


0.50 
0.17 
0.054 
0.017 
0.0054 


1.9 x 10 -4 
8.0 x 10" 6 
2.5 x 10 -7 
8.0 x 10" 9 
2.5 x 10" 10 


1.5 x 10 -30 
5.1 x 10" 31 

1.6 x 10" 31 
5.1 x 10~ 32 
1.6 x 10" 32 


1.0 
2.9 
9.2 
29 
92 



TABLE VI: Oscillaton or boson star configurations for which 
1% of the mass is lost during a Hubble time period in the 
inflationary era. For each scalar field mass m the amplitude 
parameter e, the total mass M, and the radius containing 95% 
of the mass is listed. In the last column we give the relation 
to the e = £ m ax state, which has maximal mass and minimal 
radius. 



mc 
~GcV 


h 


Mc 
MA 


T H dM 
M dt 


T H dM 
M dt 


A-0 


10 15 


0.095 


38 


0.10 


1.7 x 10 -16 


2 x 10 15 


0.047 


150 


1.5 x 10" 7 


3.3 x 10~ 16 


3 x 10 15 


0.032 


350 


5.8 x 10~ 14 


5.0 x 10~ 16 


4 x 10 15 


0.024 


620 


1.4 x 10" 20 


6.6 x 10~ 16 


5 x 10 15 


0.019 


970 


2.3 x 10" 27 


8.3 x 10 -16 



TABLE VII: Relative mass loss rates extrapolated for a Hub- 
ble time period of maximal amplitude configurations during 
the inflationary era. For each choice of scalar field mass m 
the rescaled Hubble constant h, the ratio /j c /ma, the mass 
loss rate induced by the cosmological constant, and the A = 
oscillaton mass loss rate are listed. 



likely to influence the efficiency of how the inflaton's en- 
ergy is transferred to other fields. Since oscillatons are 
concentrated energy lumps, their influence on the forma- 
tion of inhomogeneities may also be important. At the 
reheating stage the equation of state can be well approx- 
imated by a pressureless fluid, and hence the effective 
cosmological constant can be taken to be zero. For the 
lifetime of oscillatons in this era one can apply the con- 
siderations given in [35j which are valid in the A = 
asymptotically flat case. The mass of the inflaton field 
when it is near the vacuum value can be estimated as 
to = 10 13 GeV A /c 2 (see e.g. [28]). The maximal ampli- 
tude stable oscillaton which can be formed by this field, 
when e = e max ~ 0.5, has the total mass M — 0.023/c<? 
and radius r^ = 1.8 x 10 _28 m. The mass and radius of 
smaller amplitude oscillatons formed by this field can be 
easily obtained by using the fact that the total mass is 
proportional to e and the radius is proportional to 1/e. 
If the cosmological constant is zero, then a natural time 
scale is the oscillation period, T w = 2tt/uj. According to 
(f6T)|). for small e we can replace u; by to. Taking the total 
mass from (IB15|) . the part of the mass lost during an 
oscillation period is 



T u dM 

~M~dT 



Attci 

Si£ 3 



C2 

■exp I 



(140) 



This expression is independent of the scalar field mass to 
and on the units that we use for measuring time. Even for 
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the maximal amplitude oscillaton with e — 0.5 this gives 
a very tiny value, jf^jf = —2.5 x 10~ 17 . This shows 
that oscillatons formed in the reheating period perform 
a large number of oscillations before the mass loss by 
emitting classical scalar field radiation becomes appar- 
ent. However, since their oscillation period is extremely 
short, T u = 4.1 x 10~ 37 s, they lose about half of their 
mass in 8.3 x 10~ 21 s. On the other hand, since this mass 
loss depends exponentially on e, even during a period 
corresponding to the lifetime of the universe the ampli- 
tude parameter does not decrease much below e = 0.166, 
where ^l^- = -0.13. 



IX. CONCLUSIONS 

We have constructed spherically symmetric, spatially 
localized time dependent solutions in a large class of 
scalar theories coupled to Einstein's theory of gravita- 
tion and a positive cosmological constant, A. Examples 
include boson star-type objects (for a complex field) and 
oscillatons (for a real field), considered previously for van- 
ishing cosmological constant. A positive constant has 
important qualitative effects on boson stars, in that due 
to the repulsion induced by A only radiating solutions 
exist, leading to a mass loss. For sufficiently small val- 
ues of A this effect is negligible, so boson stars with al- 
most unchanged mass may be present in our Universe. 
Even when A = oscillatons radiate slowly, and the 
consequent mass loss is maximal for large amplitudes. 
There is an additional mass loss of oscillatons induced 
by a nonzero A > 0, which dominates for small ampli- 
tudes (corresponding to large sizes). The mass loss of 
an oscillaton during a Hubble time period, is completely 
determined by the ratio /i c /MA, where fi\ is the energy 
density due to A and fi c is that of the oscillaton in its 
center. Our computations are valid in spatial dimensions 
2 < D < 6, although the case D = 3 has been given spe- 
cial attention. An important cosmological application of 
our results is for " waterfall" -type oscillatons appearing 
in hybrid inflation models. 



Appendix A: Connection formulae 



In this appendix we derive the formulae (|104p connect- 
ing the amplitudes in the y < 1 and y > 1 regions. We 
follow the method described in [50]. Since the expres- 
sions for Z obtained by the WKB method are singular 
at y = 1, it is necessary to find a solution of (1751) which 
is valid in a region around this point. Let us introduce a 
new radial coordinate by 

y = l + x. (Al) 

We approximate p 2 for small h by 

p 2 w y 2 - 1 = x (x + 2) w 2x . (A2) 



Introducing a rescaled coordinate z by 

h 2/3 



751) takes the form 



d 2 Z 
dz 2 



'2V3 ' 



zZ = 0. 



(A3) 



(A4) 



The solution can be written in terms of Airy functions 



Z = a Ai z + bHi z . 
For z 3> this has the asymptotic form 



(A5) 
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^exp(--z 3/2 J +bexp(-z 3/2 j 


" V^zV4 


i for z < 


1 f _ 


-{-zf' 2 + - 
.3 4 . 




v^(-2) i/4 r 




+ b cos 


2 ( _ z) 3/2 + ^ 
.3 4 . 


}■ 






) 



(A6) 



(A7) 



We consider two matching regions surrounding y = 1, 
where \z\ 3> 1 and at the same time \x\ <C 1. Both of 
these can hold simultaneously if h is small enough. First 
we continue through the matching region with x < 0. 
Substituting y = 1 + x into (|91[) and expanding in power 
series around x = 0, we obtain 



i f v 1 3 / r \ 

W Pady=—{-2x)*+0[x*) 



(A8) 



This result can also be obtained by simply approximating 
po by v2x. In this region the influence of the si term is 
subleading, so for 3 < D < 5 dimensions we obtain both 
from dHU) and ® 



Z = 



A, 



■ exp 



(-2a;) 1 /4 
Comparing with (|A6|) . we get 

a=0 , b= 



— {-2xf 2 



(2fc)V« 



A, 



(A9) 



(A10) 



Continuing into the second matching region with x > 
0, dM} yields 



Z = 



2 y?(-z)VM ex nT' exp 



3 V 7 



+ CX P I —£ I ex P 



For small positive x (|92|) gives 



^{-z? /2 



(All) 



iL**y-wiM i+0 (**)- (A12) 
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To leading order both (|90[) and (jlOOl) can be written as 



B± 



■ exp 



(2a:) 1 / 4 
Comparing with (|A11[) gives 



±i^ 3 ' 2 



b(2h) 1 / 6 ( in 



(A13) 



(A14) 



Substituting b from (|A10[) yields the desired formulae 

mm . 



Appendix B: Mass in spherically symmetric 
asymptotically de Sitter space-times 

Since the scalar held tends to zero exponentially, 
at large distances the metric approaches the static 
Schwarzschild-de Sitter metric. Considering D + 1 di- 
mensional space-times, in Schwarzschild area coordinates 
the Schwarzschild-de Sitter metric has the form 



ds 2 



-F{f)dt 2 



1 



F(f 



-df 2 



2 dfl 



D-l ) 



where 



„D-2 



F(f) = 1 - 



-.D-2 



~H 2 f 2 



and tq is a constant related to the mass M by 



M 



{D - 1)ttt 

8-r (f ) 



n-2 



(Bl) 



(B2) 



(B3) 



In general spherically symmetric space-times there is 
a naturally defined radius function r, dehned in terms 
of the area of the symmetry spheres. In Schwarzschild 
coordinates f = f, while in the isotropic coordinates (|19[) 
f = ryB. In terms of this radius function one can dehne 
the Misner-Sharp energy (or local mass) function rh [5l|, 
|52j , which can be defined for arbitrary dimensions by 



(D-I)tt^ 

**r (?) 



lD-2 



(1 



<A,„r 6 



# J r 



2~2\ 



(B4) 



For the Schwarzschild-de Sitter metric rh = M . At in- 
finity the mass function rh agrees with the Abbott-Deser 
mass [53l . I54J , which is a Killing energy defined in static 
asymptotically de Sitter space-times. 

The Misner-Sharp energy rh also agrees with a con- 
served energy that can be defined using the Kodama vec- 
tor. The Kodama vector [55|, [56J is defined by 



K a 



e ab f,b , 



(B5) 



where e a b is the volume form in the (t, r) plane. Choosing 
the orientation such that e rt — \AB makes K a future 
pointing, with nonvanishing components 



K< 



CAB 



K r 



>AB 



(B6) 



It can be checked that, in general, the Kodama vector 
is divergence free, K a a = 0. Since contracting with the 
Einstein tensor, G K a; i, = 0, the current 



J a = T ab K 



(B7) 



is also divergence free, J a a = 0, it defines a conserved 
charge. Integrating on a constant t hypersurface with 
a future oriented unit normal vector n a , the conserved 
charge is 



E 



2ir~ 






r(#)./o 



f D - 1 y/Bn a J a dr (B8) 



r(#)V A 



(T tt r, r - T t rf } t) dr . 



It can be checked that the derivative of the mass function 
rh can be expressed in terms of the current J a , 



2tt 2 f 



#££>-! 



m.a = -- 



r(#) 



-tab J 



For the radial derivative follows that 

(Tur^r - T tr r,t) , 



r(f)A 



(B9) 



(BIO) 



which, comparing with (IB8|) . gives E = rh. 
The time derivative of the mass function is 

2^f D - 1 , 
m, t = wrm p \ T rtr, r - T rr r, t ) ■ 



r(f)B 



(Bll) 



This equation describes the mass loss caused by the out- 
ward energy current of the massive scalar field. 

For large radii the function rh tends to the total mass 
M of the oscillon or boson star. Substituting the small- 
amplitude expansions (|34D and (|35l) into the definition 
(|B4[) of the mass function m, to leading order in e we 
obtain 



_ 4 _ D 7TT(D-i) p g-i dB 2 
8Trm D - 2 r(%) dp 



(B12) 



We have seen in Subsections IIII B land MI CI that for both 
the boson star and oscillaton case B 2 = A 2 /(2 — D) and 
A% = — 1 — s, giving 



-4-D 



m 



Tr-p-l)^- 1 ds 



8ir(D - 2)m D - 2 T (f ) dp' 
Since the asymptotic behavior of s is given by ()58[) , 



(B13) 



87rm^- 2 r(f) Sl ' 
For 3 + 1 dimensional space-time 

M(r>=$\ = £-z — ■ 
Am 



(B14) 



(B15) 



For small h the numerical value of the constant si can 
be approximated by the h — value given in Table [TTJ 
Measuring the mc 2 belonging to the scalar field in elec- 
tron volts and the mass of the boson star or oscillaton in 
kilograms, 



eV 



M {D=3) = 4.657 x lO^kg — ^s 



mc* 



(B16) 



Appendix C: Proper mass and core radius 
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D = 3 


D = 4 


73 = 5 


P50 


2.240 


1.778 


1.317 


P'JO 


3.900 


3.013 


2.284 


P95 


4.471 


3.455 


2.652 


P99 


5.675 


4.410 


3.478 


P99.9 


7.239 


5.692 


4.634 



TABLE VIII: The radius inside which given percentage of the 
mass-energy is contained for various spatial dimensions. 



The proper mass inside a sphere of radius r is defined 
by the D dimensional volume integral of the energy den- 
sity /i, 



M p (r) 



2tt* 

r(f) 



drpB 



D/2 r D- 



(Cl) 



Substituting the leading order expression (|64[) for p, and 
the leading order value B = 1 from (j35|) . using the 
rescaled radial coordinate p = emr, 



M p (p) 



-4-D 



7TT(D-1) 



8n(D - 2)m D - 2 T (§) J 



dpp u - l S 2 . (C2) 



The integral can be performed using the SN equation 
(El, 



M p (p) 



e 4 - u Tr^(D-l)p 
iir(D - 2)m D - 2 r 



D-l 



. 2 / 



ds 

dp 



(C3) 



To leading order in e the calculated proper mass M p 
agrees with the mass function m given by (JB13I) . How- 
ever at higher order the proper mass is expected to be 
larger by an amount corresponding to the binding energy. 
Oscillatons and boson stars do not have a definite outer 
surface. A natural definition for their size is to take the 
radius r n inside which n percentage of the mass-energy 
can be found. It is usual to take, for example, n = 95 
or 99.9. The mass-energy inside a given radius r can 
be defined either by the integral (IC1|) or by taking the 
local mass function rh in (|B4I) . To leading order in e 
both definitions give (|C3j) . The rescaled radius p n can 
be defined by 



Mp(p n ) _ _n_ 
Mp(oo) ~ 100 ' 



(C4) 



The numerical values of p n for various n in D = 3,4,5 
dimensions are listed in Table IVIIII The physical radius 

is 



We note that for complex field boson stars the current 



J 



2 9 



$*$,-$*$ 



(C7) 



can be used to define a conserved quantity. The con- 
served particle number N can be defined by integrating 
the time component j° for a spacelike slice. To leading 
order in e the function j° is proportional to the density p 
given in (IM1) . Hence, to this order, the use of the particle 
number instead of the mass-energy yields the same result 
for the radius of boson stars. 



Appendix D: Present value of the cosmological 
constant 

For the present observational value of the Hubble con- 
stant we take 

H a = 70.4—^- = 2.28 x 10~ 18 - (Dl) 

sMpc s y ' 

according to the combined 7- year WMAP data [53|. To 
this corresponds a critical density 



l^cr 



8irG 



= 9.31 x lO- 24 -^- 



(D2) 



The fraction of this belonging to the dark energy pro- 
vided by the cosmological constant is given by 51a = 
0.728, according to WMAP results [57|. Consequently, 
the energy density belonging to the cosmological constant 

is 



p A = p cr fl A = 6.78 x 10~ 24 -^- 



(D3) 



According to (|31[) . this corresponds to a cosmological 
constant of 



8irGu\ n 1 

A = J-± = 1.26 x 10~ 52 ^ 



(D4) 



em 



(C5) 



In ordinary units, measuring the scalar field mc 2 in elec- 
tron volts and r„ in meters (Roman m), 



1.97xKr 7 m^^C 



e mr 



(C6) 



In Planck units, since the Planck length is lp — 
1.616252 x 10~ 35 m, 



A = 3.30x 10" 122 . 



(D5) 



Since we are interested in the effect of the cosmologi- 
cal constant, we neglect the ordinary and dark matter 



18 



content of the universe. Then the space-time can be de- 
scribed by a de Sitter metric with Hubble constant H 
given by (|26[) . In Planck, ordinary and natural units 

H = 1.05 x 1(T 61 = 1.95 x lCT 18 i = 1.28 x l(T 33 eV. 

(D6) 
The Hubble time is then 



1 



T H = — =5.14x l(r's = 16.3Gyr. (D7) 

H 

Even though the ordinary and dark matter content have 
been neglected, this value is still not too far from the 



age of the universe, which according to WMAP data is 
13.75Gyr assuming the ACDM model [57J . 
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